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D : Abstract 

Effects of a periodic driving field on Landau-Zener processes are studied using a nonlinear two- 
C^ ■ mode model that describes the mean-field dynamics of a many-body system. A variety of different 

dynamical phenomena in different parameter regimes of the driving field are discussed and ana- 
^ ' lyzed. These include shifted, weakened, or enhanced phase dependence of nonlinear Landau-Zener 

o ■ 

, ^, ' processes, nonlinearity-enhanced population transfer in the adiabatic limit, and Hamiltonian chaos 

fSJ ■ on the mean field level. The emphasis of this work is placed on how the impact of a periodic driv- 

>, 
OO ing field on Landau-Zener processes with self-interaction differs from those without self-interaction. 

m 

•^ ' Aside from gaining understandings of driven nonlinear Landau-Zener processes, our findings can be 

T^ij- ■ used to gauge the strength of nonlinearity and for efficient manipulation of the mean-field dynamics 

o. 

OO , of many-body systems. 

O 
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I. INTRODUCTION 

Landau-Zener (LZ) transition [1] forms a fundamental dynamical process relevant to 
a variety of contexts, such as quantum control and quantum information [j, y, 14, y, |6|, 



B 



quantum dots [7[, molecular clusters [8j, to name a few. If the external bias involved in a 
LZ process is tuned adiabatically, then the LZ tunneling probability goes to zero, yielding a 
robust scenario for realizing a complete quantum population inversion. If the external bias 
is varied non-adiabatically, then the LZ tunneling offers a diagnosis tool for understanding 
the quantum dynamics. Indeed, the paramount importance of LZ processes has motivated 



a large body of theoretical and experimental work [9|, [lOj. 

As a very recent development, Wubs et al. [11] studied numerically and analytically how 
a LZ process might be further manipulated by considering a periodic driving field. Within 



the rotating- wave-approximation (RWA), Wubs et al. [ll| showed that a LZ process in 
the presence of a driving field possesses two intriguing features. First, the driving field can 
induce interesting quantum interference effects between two well-separated LZ sub-processes, 
with the final quantum transition probability dependent upon the phase of the driving field. 
Second, contrary to a standard LZ problem, the population inversion probability approaches 
zero (instead of unity) if the bias is tuned adiabatically. These findings may find novel 
applications of LZ processes, especially in identifying an unknown phase of a driving field 
or in studying decoherence time scales of a system coupled with an environment. 

Generalizing LZ processes for single quantum systems to those associated with the mean- 
field dynamics of interacting many-body systems, one obtains nonlinear LZ (NLZ) processes. 
Due to their implications for quantum control of the dynamics of Bose-Einstein condensates 



(BEG) 12|, NLZ processes have also attracted considerable interests. Note that the mean- 
field dynamics of an interacting many-body system is necessarily nonlinear, and as a result 
simple physical intuitions based on linear LZ processes may become invalid in nonlinear 
cases. For example, the adiabatic following of a quantum state with an external bias may 
)reak down in an NLZ process, even when the external bias varies at an infinitely slow rate 
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Extending the work by Wubs et al. [ll| from the linear regime to the nonlinear regime, 
our interest here is in how a periodic driving field might affect NLZ processes. It is hoped 
that a complete picture of driven NLZ processes will help make the best use of a driving field 



in manipulating the mean-field dynamics of many-body systems. To reach that long-term 
goal, it is necessary to first examine how the impact of a periodic driving field on NLZ 
processes differs from that on conventional LZ processes. 

In our previous work [16], NLZ processes subject to a high-frequency and large-amplitude 
driving field were studied. There, using a high-frequency approximation, we showed that 
the driven NLZ dynamics can be effectively described by a stationary Hamiltonian. The 
spectrum of the effective Hamiltonian can have new degenerate eigenstates and hence display 
new topological structures that are absent in the non-driven cases. The results offer a simple 
approach to the complete suppression of NLZ tunneling, even when the tunneling is doomed 
to happen in the absence of a driving field. 

This work continues to study driven NLZ dynamics, but with the parameter regimes 
different from that in Ref. 16|]. We place our emphasis on comparisons between driven 
(linear) LZ processes and driven NLZ processes, thus shedding more light on the nonlinear 
and mean-field nature of NLZ dynamics. For example, we shall show that nonlinearity can 
induce significant quantum population transfer that is completely absent in linear cases. We 
will also show that the extent of population transfer in driven NLZ processes can be, as 
compared with driven LZ processes, either much more or much less sensitive to the phase 
of the driving field. It is hoped that our detailed results below will motivate more efforts 
towards better understanding and better control of many-body systems. 

This paper is organized as follows. In Sec. II we briefly discuss our two-mode model 
of NLZ processes in a periodic driving fleld. In Sec. Ill we discuss the high-frequency 
approximation in treating NLZ processes in a high-frequency and large-amplitude driving 
fleld. In Sec. IV we consider a different regime of fleld parameters, where the RWA can be 
applied. In Sec. V we study another parameter regime, where neither the high-frequency 
approximation nor the RWA is valid. Section V concludes this work. 



II. MODEL OF NONLINEAR LANDAU-ZENER PROCESSES IN AN "OFF- 
DIAGONAL" DRIVING FIELD 



Motivated by the model of linear driven LZ processes considered in Ref. ll|], we consider 
a two-mode model of driven NLZ processes as follows: 



Here, 



2 \ Ac + Aocos(c^t + /3) -7-c(|6|2-|a|2) 
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7 = at 
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2l|). the driving field 



denotes a time-dependent bias between two modes of interest and is being varied at a rate 
a; |ap and |6p represent occupation probabilities on the two modes, with the normalization 
condition taken as |ap + |6p = 1. The terms containing c characterize the nonlinear self- 
interaction of a BEC under the mean-field treatment, with the value of c proportional to 
the number of bosons and the s-wave scattering length |l2l . Il3l . Il4l . [l5| ; Ac represents the 
static coupling between the two modes; and Aq is the amplitude of an external driving, with 
the frequency u and the phase parameter 3. Unlike previous models for driven LZ/NLZ 
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processes pJl, Il7|, Il8| or for coherent destruction of tunneling [9|, M (fo r a recent beautiful 
experimental validation of the latter phenomenon see in Refs. 
here appears only in the off-diagonal terms of the above Hamiltonian and is hence an "off- 
diagonal" driving field. That is, the external field directly modulates the coupling between 
the two modes, rather than modulating their energy bias. The initial state of an NLZ process 
is always taken as a = 1 and 6 = 0, and the occupation probability |6p in the end, denoted 
by |6(oo)p, is the final quantum transition probability (or the final population transferred 
between the two modes). Note also that all the variables here should be understood as 
scaled dimensionless variables with h = 1. Throughout we use Aq to scale. Then, Aq = 1, 
uj is in units of Ao/h, a is in units of A^/h, and c is in units of Aq. 

Within the BEC context, the above Hamiltonian may be experimentally realized in several 



ways. For example, one may consider a BEC in a double-well potential 22 



23 



24 



25|, with 



the barrier height periodically modulated, or by a BEC in an optical lattice occupying two 
bands H. with the well-depth of the opfeal latfce per.od.cally modulMed. One may also 

regard the two modes as two internal states of a BEC, such as ^^Rb [27j, with the energy 



bias 7 effectively realized by the detuning of a coupling field from the resonance and the 
off-diagonal modulation achieved by modulating the intensity of the coupling field. Though 
we only refer to the BEC context below, it should be noted that this two-mode model of 



NLZ might be relevant to physics of Josephson junctions as well as nonlinear optics [/ 

In the special case of Aq = 0, the above Hamiltonian reduces to the well-known model 
of non-driven NLZ processes [ij]. Therein the eigen-spectrum diagram as a function of 7 is 
known to display a loop structure at the tip of the lower (upper) level for c > Ac (c < —Ac)- 
Such a loop structure, absent in linear systems, directly leads to a nonzero LZ transition 
probability even when 7 changes adiabatically. This makes it interesting to examine what 
happens if an NLZ process is subject to a periodic driving field as introduced above. 

Without loss of generality we will restrict ourselves to the c > cases (cases with c < can 
be mapped to those with c > 0). This c > assumption requires, for example, an attractive 
interaction for bosons in a double-well potential or a repulsive interaction for bosons in two 



energy bands of an optical lattice. In current non-modulated double-we 
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]C experiments, 



25|. For the latter 



typical values of |c/Ac| for 10 bosons range from ~ 10 to ~ 10 

two-band realization of our two-mode model, the value of c (scaled by An) can be easily 

tuned by controlling the maximal depth of an optical lattice potential 14 1. 

To examine the driven NLZ dynamics (i.e., Aq 7^ 0), we consider below three different 
parameter regimes. In the first high-frequency and large- amplitude regime , a high-frequency 
approximation can be used to obtain an effective Hamiltonian. The second regime is the 
RWA regime, where the driven dynamics is again understood in terms of some effective 
Hamiltonians. In the third regime, both the driving frequency u and the amplitude Aq are 
comparable to, or smaller than, the nonlinear parameter c. Our task below is to examine 
the driven NLZ dynamics in these three regimes. 

III. HIGH-FREQUENCY AND LARGE- AMPLITUDE DRIVING FIELD 

The NLZ dynamics under the condition a;,Ao ^ 70, c, Ac, 70 (where 70 is the initial 



value of 7) was first studied by us in Ref. 16j using a high-frequency approximation. For 



completeness and for the sake of comparison with other regimes, here we briefly revisit this 



regime. To that end we first introduce another pair of wave function parameters (a', b') 

, _ Q + ^ i4iasinM+/3) a--b_ ^^^i^(^^^^p) 
2 2 

Substituting these two relations into the Hamiltonian in Eq. ([T]), one finds the equations of 
motion for (a', b'): 

,^ = l[^.cos(e) + c-cos2(e)(|6f -|af) 
+ic ■ sm{9) cos{9){a*b' - ab'*)]a' 
+-[Ac - ii ■ sin(^) + c ■ sin2(^)(a'*6' 



-a'6'*) + %c ■ sin(^) cos(^)(|a'|^ - |6f )]6' 

2^ = i[Ac; + 27-sin(e)-c-sin2(^)(a'*6' 

-a'r) - ic ■ sin(^) cos(^)(|af - |6f )]a 
1 



f -|6f)]a' 

2//i\/lr,/|2 |„/|2\ 



+-[-7-cos(^)-c-cos'(e)(|6' 

-ic ■ sin(^) cos(^)(a'*6' - a'6'*)]6', (4) 

where ^ = ^ sin (cut + /5). 

For sufficiently large a;, the oscillation in ^ is much faster than all other time scales of the 
system. As such, during each period of T = 1i\ juj, the change in o! and b' is negligible and 
considering their averages over T will greatly reduce the equations of motion. Specifically, 
upon a time averaging the odd functions of t in Eq. (j4]) vanish and the effective equations 
of motion become 

• r ' I /i!,'i2 I n2M / 

^"dF " T^ +cz(|fe| -|a| )]a 

+ ^[Ac + cy(a'*6'-a'r)]6' 

1% = hAc - cyia'*b' - a'b'*)]a' 
at 2 

+ l[-i-cz{\bf-\af)]b', (5) 

where 7' = 7(cos(e))T = 7^o(Ao/cu), cz = c{cos^{e))T = c[l + Jo{2Ao/uj)]/2, cy = 
c{sm'^{9))T = c[l — Jo(2Ao/ci;)]/2, and Jq is the zeroth order Bessel function of the first 
kind. Apparently, these newly defined parameters reflect the action of the high-frequency 
driving field. 
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FIG. 1: Eigen-energy structures of the effective Hamiltonian in Eq. ([6]) that describes NLZ pro- 
cesses in a high-frenquency and large-amphtude driving field. The effective Hamiltonian is obtained 
with a high-frequency approximation. The discrete points represent the time-evolving expectation 
value of the mean-field energy when 7 increases in one adiabatic process. Along with other parame- 
ters, the values of the nonlinear parameter c and the driving frequency u are indicated on the three 
panels. As mentioned in the text, all variables are appropriately scaled and hence dimensionless. 



Based on Eq. ([5]), we can define the following effective Hamiltonian: 



H. 



eS 



7 + cz(|&p-|an Ac + CY{a*b-ab 
Ac - CY{a*b - ab*) -7 - cz(|6p - 




(6) 



where we have replaced a' by a, b' by 6, and so on. Comparing this effective Hamiltonian 
with the original one in Eq. ([T]) for Aq = 0, it is seen that the nonlinear parameter Cz can 
be regarded as a re-scaled parameter c, and the new nonlinear terms containing cy arise 
as a surprise. In addition, the ratio of the two nonlinear parameters cz and cy is given by 
[1 + Jo(2Ao/u;)]/[l — Jo(2Ao/u;)], which is easily adjustable by tuning the ratio Aq/oj. 

In terms of the mean-field eigen-energies of iifefr as a function 7, three typical level struc- 
tures of i^eflf are depicted in Fig. 1. The structures represented by the solid lines in Fig. 1 are 
also typical in non-driven NLZ models [14], but the dashed lines (each of them is associated 
with double degenerate eigenstates) in Fig. 1 are entirely induced by the high-frequency 
driving field. For an NLZ process with 7 increasing at a very small rate of a = 0.002 (hence 
an adiabatic NLZ process), the time-evolving expectation values of the mean- field energy 



associated with Hes are shown by the discrete points in Fig. 1. Clearly, for the two cases in 
Fig. 1(a) and 1(c), the evolution of the system perfectly follows the new eigenstates (dashed 
lines). The driving field can hence dramatically affect the dynamics of adiabatic following. 
Indeed, in the case of Fig. 1(a), if there were no additional levels induced by the driving 
field, then the loop structure there will necessarily cause the adiabatic following to break 
down (this is analogous to what is observed in non-driven NLZ models [ij] and will be also 
seen below). For other important findings in this high-frequency large-amplitude regime, 



please see Ref. [161]. 



The high-frequency approximation used here is also applicable to high-frequency and 
small-amplitude fields. However, if Ao/cj << 1, then Jo(2Ao/ct;) ~ 1, c^ ~ c, cy ~ 0, and 
the effective Hamiltonian in Eq. (6) will essentially reduce to the Hamitonian of Eq. (1) 
with Ao = 0. That is, in cases of Aq/uj << 1, the present high-frequency approximation 
based on the (a', b') representation does not directly give useful insights into the dynamics. 
Fortunately, in these cases the RWA treatment becomes valid and more advantageous to 
use. This is discussed in detail in the next section. 

IV. HIGH-FREQUENCY AND SMALL- AMPLITUDE DRIVING FIELD 

A. Rotating- Wave Approximation 

In the high-frequency and small-amplitude regime, we adopt the RWA to understand the 
effects of a driving field on NLZ processes. The off-diagonal coupling term in Eq. (1) is now 
understood as a superposition of three terms, 

A = — [exp{-iujt - i(3) + exp{iujt + i(3)] + Ac, (7) 

where the first two terms represent the circular modulation with frequency uj and —u, 
and the last constant term can be regarded as a circular modulation with zero frequency. 
Consistent with the high-frequency assumption, the effects of the three circular modulation 
terms can be analyzed separately. To illustrate this, we depict in Fig. 2 a typical numerical 
result for this parameter regime. Evidently, the transitions primarily occur in three small 
time windows around 7 ~ —uj, 0, u, each of them is associated with one of the three circular 
modulation terms in Eq. (jTj). 
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FIG. 2: Transition probability |6(t)p as a function of at, with the phase parameter (3 = 0. The 
driving field is in the RWA regime. Three well-separated NLZ transitions can be seen here. 



For the sake of a direct comparison between driven NLZ processes and the driven LZ 
processes studied in Ref. [ll|, below we set A,^ = 0. We then only need to consider two 
circular modulation terms. The effects of the first circular modulation term -^ exp(— ■jcut— i/?) 
become obvious in a rotating frame. We hence introduce another pair of wavefunction 
parameters (a, b). 



aexp 
6exp 



--t{ujt + f3) 



-Aivt + p) 



(8) 



Using this representation and neglecting the counter-rotating terms, one obtains 



.d_ 




H+{t) 




(9) 



where the RWA effective Hamiltonian H^{t) is given by 






Ao 
2 



lap) 



(10) 



with 7+ = 7 + cu. This Hamiltonian is essentially identical with that for a non-driven NLZ 
model, with the bias parameter 7 shifted by 00 and with the effective static coupling term 



determined by the amplitude of the periodic driving field. This explains the observation 
from Fig. 2 that one main transition window is roughly at 7 ~ — cu. Because |6p = \b\^, 
this RWA treatment also indicates that for a fixed initial state, the quantum population 
at the end of such a transition window should not depend on (3. Similar treatments apply 
to the other circulation modulation term ^exp(iu;t + iP), yielding the following effective 
Hamiltonian 

^\ t -1. - c(\bf - \af) ) ' 

with a' = aexp [^^(cut + P)~\ and b' = 6exp [—ji{ujt + /3)]. Here the effective bias is 7_ = 
'-f — uj, and as a result the main transition window will be close to 7 ~ tu. 

Because the RWA treatment here is seen to be independent of the nonlinear parameter 
c, it applies as well to those driven (linear) LZ processes considered in Ref. |11|]. Hence, for 
either NLZ or LZ processes, a periodic driving field is expected to shift the main quantum 
transition window at 7 ~ to two new windows close to 7 ~ ^u (how the exact locations 
of the transition windows depend on c will be discussed later). 

B. Phase-dependence w^eakened by strong nonlinearity 

The above RWA treatment makes it clear that the NLZ dynamics in this regime can 
be described by effective Hamiltonians analogous to that for non-driven cases. However, 
because the final state of the transition window around 7 ~ —u becomes the initial state 
of the second transition window around 7 ~ cj, the driven dynamics should be much richer 
than in non-driven cases. Indeed, in the linear case with c = 0, the analytical result from 



Ref. [HI] predicts that the final transition probability is given by 



|6(oo) |2 = 4 exp(-7rA^/8«) [l - exp(-7r Ag/8a)] cos^ (3, (12) 

where (3 is the phase of the driving field defined above. Because the transition probability 
of the first LZ sub-process does not depend upon (3, the factor cos^(/5) in Eq. flT2l) refiects 
how the net result associated with the second transition window is affected by the first one. 
Alternatively, this cos^(/?) factor can be interpreted as a quantum interference between two 



LZ sub-processes at two time windows HI]. This intriguing phase dependence might be 
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FIG. 3: Left panel: Transition probability |6p as a function of at, with a = 0.56, lo = 37.5, and 
different values of the nonlinear parameter c and the phase parameter (3. Right panel: The final 
NLZ transition probability |5(oo)p (solid lines) as a function of f3. Dashed lines represent the 
theoretical linear result of Eq. ()12p that displays a cos^(/3) dependence. 



useful in determining an unknown value of /3 11|. 



How does the nonlinear self-interaction of a BEC affect this cos^(/3)-dependence of the 
final transition probability? Unable to obtain an analytical solution, we choose to examine 
this issue computationally. Typical results for a rather rapid ramping of 7 are shown in Fig. 
3. The ramping rate is chosen to be a = 0.56 such that the linear result of Eq. (TT2l) yields 
the maximal transition probability. Panels (a), (b) and (c) of Fig. 3 depict the time-evolving 
occupation probability |6p as a function of at, with the initial condition a = l,b = 0. In 
each of the three panels, results for a number of f3 values are plotted collectively. For all 
the shown cases two main quantum transition windows can be clearly seen, consistent with 
our expectation based on the RWA analysis. Comparing Fig. 3(a), 3(b) with Fig. 3(c), it is 
seen that as the nonlinear strength c increases, the final transition probabilities for different 
values of (3 tend to converge to a common value, i.e., they become less and less sensitive 
to p. To be more quantitative, we extract the /3-dependence from Figs. 3(a)-3(c) and then 
compare it (right panels in Fig. 3) with the linear result of Eq. ( IT2l) . For the case of c = 0.25 
(right upper panel of Fig. 3), the /3-dependence of the final transition probability |6(oo)p 
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exhibits a remarkable phase shift as compared with the hnear result (dashed line). This 
phase-shift effect might be useful in determining a small but unknown value of c. For cases 
with sufficiently strong self-interaction (right bottom panel of Fig. 3), |6(oo)p changes little 
as (3 varies, constituting a sharp contrast to the cos^(/5) dependence in linear cases. 

The weakened /5-dependence due to strong nonlinearity requires a qualitative explanation. 
Consider the state after the first transition window, i.e., the initial state for the second NLZ 
sub-process. As (3 varies, the occupation populations associated with this intermediate state 
can be hardly changed (as is evident from the above RWA treatment), but the relative phase 
between the amplitudes a' and b' will be affected. With this understanding, we explain below 
the observed weakened /^-dependence of |6(oo)p in terms of the increasing insensitivity of 
the second NLZ sub-process on the quantum phase of its initial state. 

Let us map exactly the mean-field dynamics of the effective Hamiltonian in Eq. flTTj) to 
that of a classical Hamiltonian system, with the canonical variables S = |6'P — la'P and 
Il8l |. with (J) the relative phase between the amplitudes a' and b'. The mapped 



y, y, L 



classical Hamiltonian, denoted by He, is then given by |l3| : 

1 



He 

2 



j_S-^S' + ^Vl^S^cos, 



(13) 



With this mapping, the phase of the initial state for the second NLZ sub-process apparently 
becomes the initial value of the coordinate 0. Further, the final transition probability |6(cxd)P 
is given by [5'(oo) + l]/2. 

Figure 4 compares the phase space structure of He for 7_ = 7 — cj = and that for 
7_ = 100. Consider first the upper two panels for c = 0. In the regime of large 7_ [Fig. 
4(a)], the phase space curves become straight lines perpendicular to the 5* axis, with their 
S values determining |6(oo)p. In the regime around 7_ = [Fig. 4(c)], the phase space 
invariant curves rotate around elliptic fixed points, with their oscillation amplitude in S 
sensitively dependent upon 0. Clearly then, if 7_ is scanned rather rapidly, the range of 
S values for 7_ = 0, which depends on the initial value of 0, will determine the range 
of S and hence the possible values of |6(oo)p. This qualitatively explains the sensitive 
/3-dependence in the linear case. By contrast, as shown in Fig. 4(b) and Fig. 4(d), the 
phase space structure for large c is largely independent of 7_, with most of the phase space 
curves being parallel fines perpendicular to the S'-axis. The sensitive dependence of S{oo) 
on 0, and hence the sensitive dependence of |6(oo)^| on j3, must be weakened for large c. 
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FIG. 4: Phase space structure associated with the mapped classical Hamtonian in Eq. (J13p . in 
terms of the dimensionless variables S and (j) defined in the text. For large c, most phase space 
invariant curves will be, independent of the value of 7_ = ^ — lo, almost straight lines perpendicular 
to the S'-axis. 

Indeed, in the large c limit, one may only keep the —cS'^/2 term of He- He then becomes 
the Hamiltonian for a "free-particle" system, with a trivial phase space structure filled with 
straight lines perpendicular to the S'-axis. As such, the /3-dependence of |6(oo)^| in the large 
c limit will disappear altogether. This completes our qualitative explanation of the weakened 
/9-dependence due to increasing nonlinearity. 

Finally, we discuss another interesting observation made from Fig. 3. Comparing Fig. 
3(a), 3(b) and 3(c), one notes that the first transition window is delayed more and more as 
the nonlinear strength c increases. The delay of the second transition window as a function 
of c is however less apparent. To rationalize this observation we first note that in the cases 
with negligible c, the standard picture based on linear LZ processes applies and the two 
transition windows should be centered precisely at 7 ~ —uj and 7 ~ +uj. In cases with large 
c, the effective Hamiltonians H^ and if_ obtained in our RWA treatment are expected to 
display loop structures, in the same manner as in non-driven NLZ models [ij]. One typical 
example with c = 17.5 and cj = 50 is shown in Fig. 5, where both i7_|_ and H^ display a 
loop structure in the lower branch of the energy levels. Due to the loop structure of if+, 
the increase in 7 from a very negative value will not cause any significant transition to 
the upper branch at 7 ~ —u. Rather, significant population transfer occurs only until the 
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FIG. 5: Schematic energy level structures of the RWA static Hamiltonians H^ in Eq. (jlOp and H- 
in Eq. (jll|) . c = 17.5 and uj = 50. All variables are appropriately scaled and hence dimensionless. 

system reaches the right "edge" of the loop, thus postponing the first transition window. In 
addition, because the loop structure is bigger for larger c, the first transition window will 
be postponed more as c increases. 

Interestingly, after the first NLZ sub-process, the system will in general occupy both 
upper and lower levels considerably. Then, as 7 continues to increase in the positive regime, 
those populations on the upper branch will make significant transitions at 7 ~ cj, whereas 
those on the lower branch will delay their transitions again. The overall result is that the 
precise location of the second transition window becomes less definitive and more dependent 
on the details of the dynamics. With this recognition we schematically place the second 
transition window in Fig. 5 between 7 ~ cj and the right edge of the right loop structure. 

C. Nonlinearity-induced Transition Probabilities 

Equation ( IT2l) contains another intriguing and stimulating result for driven linear LZ 
processes. That is, in the adiabatic limit a —>■ 0, the final transition probability |6(oo)p 
becomes zero for any /5. The associated physical picture is as follows. After the first 
transition window of a linear driven LZ process, the adiabatic population transfer is already 
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FIG. 6: Transition probability |6(f)p as a function of at, for a small value of a = 0.006 that 
simulates the dynamics in the adiabatic limit. Panel (a) and panel (b) are for two different values 
of c. In each case the results for a number of /3 values are plotted collectively. The large final 
values of |6p are induced entirely by the nonlinear nature of the system. 

complete. But then the system approaches the second transition window, which transfers 
all the population back to the initial mode. In other words, for linear driven LZ processes, 
the transitions associated with the two time windows 7 ~ ^uj exactly cancel each other in 
the adiabatic limit. Such a strong result for linear systems motivates us to examine driven 
NLZ processes in the same adiabatic limit. 

For two values of the nonlinear parameter c. Fig. 6 shows the time-evolving occupation 
probability |6p, using a sufficiently small value of a to simulate the dynamics in the adiabatic 
limit. In each case, results with different values of f3 are plotted collectively. In all these 
cases, it is seen that the quantum transition probabilities at the end of the first transition 
window are far from unity. Because in linear cases the adiabatic population transfer after 
the first transition window should be always unity, the incomplete population transfer shown 
in Fig. 6 after the first transition window constitutes direct evidence of considerable NLZ 
tunneling in the adiabatic limit. Analogous NLZ tunneling in the second time window can 
also be expected. As to the final transition probability |6(oo)p (final values of |6p), they 
are seen to be insensitive to [3. This phase insensitivity is similar to previous non-adiabatic 
cases with considerable c values. But more noteworthy, |6(oo)p is ~ 0.35 for c = 17.5 in 
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FIG. 7: The final transition probability |6(oo)p (solid line) as a function of the nonlinear strength 
c, for a high-frequency and small-amplitude driving field. The dashed line shows the transition 
probability right after the first transition window. 

Fig. 6(a) and as large as ~ 0.5 for c = 2.12 in Fig. 6(b). Considering that |6(oo)p should 
be precisely zero if c = 0, the large transition probabilities |6(oo)p here are induced entirely 
by the nonlinear nature of the system. In addition, comparing Fig. 6(a) and Fig. 6(b), 
one also sees that the larger c is, the more delayed the first transition window will be. This 
observation is analogous to what is found in Fig. 3. 

To further motivate interests in the nonlinearity-induced transition probabilities in driven 
NLZ processes, we calculate |6(oo)p in many adiabatic cases that cover a wide range of c 
values. A typical result, obtained with u = 40 and a = 0.003, is depicted in Fig. 7. The 
non-monotonic dependence of |6(oo)p upon c is apparent. In particular, for any c < 0.5, the 
transition probability right after the first transition window is close to unity and becomes 
essentially zero after the second transition window. This behavior is the same as in linear 
driven LZ processes in the adiabatic limit. Hence, for the RWA regime studied here, the 
threshold value for observing nonlinearity-induced population transfer at the end of the 
adiabatic process is around c ~ 0.5. Figure 7 also indicates that for very large values 
of c, |6(oo)|^ decays back to zero again. For intermediate values of c, the final transition 
probability can be significant. Indeed, the peak value of |6(oo)p as a function of c is found 
to be as large as ~ 0.95. 
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Noticing that the c value here is scaled by Aq, it is clear that by tuning the amplitude 
of the driving field one may realize the peak value of |6(oo)p observed in Fig. 7, without 
actually tuning the natural scattering length or the number of bosons of the BEC. Hence a 
nearly complete population inversion is made possible by a periodic driving field in the RWA 
regime. This finding offers a potentially useful scenario for improving the adiabatic popu- 
lation transfer in NLZ processes. Interestingly, in the high-frequency and large-amplitude 
regime discussed in Sec. Ill, the adiabatic population transfer is improved by suppressing 
the NLZ tunneling; whereas the nearly complete population transfer here is based on the 
overall effect of two NLZ tunneling processes. 

Combining (i) the shifted /9-dependence of |6(oo)p for small c, (ii) the weakened (3- 
dependence of |6(oo)p for large c, (iii) the delay of the first transition window as a function 
of c , and (iv) the nonlinearity-induced transition probabilities in the adiabatic limit, one 
may also design strict tests to check the validity of a two-mode description of a driven 
many-body quantum system. 

V. LOW-FREQUENCY DRIVING FIELDS 

If the driving frequency oo is comparable to, or smaller than, c and Aq, we call it a low- 
frequency case. In this regime the previous high-frequency approximation or RWA does not 
apply. As it turns out, complications arise here due to the emergence of Hamiltonian chaos 
on the mean-field level. 

Using the same approach in Sec. IV-B, let us first map the driven NLZ dynamics described 
by Eq. (1) (with A^ = 0) to that of a classical Hamiltonian, 

H^it) = - \--fS - -S^ + Ao cos(cc;t + P)VT^^cos </)] , (14) 

where S is given by (2|6p — 1), and is now the relative phase between the amplitudes 
a and b. Because u is small, it is impossible to divide the driven NLZ process into two 
sub-processes. But this also implies that, roughly speaking, the phase space structure at 
7 ~ will be most important in understanding the mean-field dynamics. As such, we show 
in Fig. 8 the typical behavior of the Poincare surface of section of Hc{t) with 7 = 0, a; = 1, 
/9 = 0, and with the nonlinearity parameter c ranging from to 250. It is seen that for 
intermediate values of c, the phase space is completely chaotic, whereas for small or large 
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FIG. 8: Phase space structure associated with the classical effective Hamiltonian in Eq. (14) with 
/3 = 0, in terms of s and (p defined in the text, for various values of c. The driving frequency is given 
by cj = 1. As mentioned in the text, all variables are appropriately scaled and hence dimensionless. 

values of c, the dynamics is primarily regular. 

In Fig. 9 we plot |6(oo)p as a function of j3, with the values of uj same as in Fig. 8. 
Consider first Fig. 9(a) for the c = case. The /9-dependence is seen to be rather weak, in 
contrast to the cos^(/3) result in the RWA regime. In addition, |6(oo)p is small for the entire 
range of f3. Next, in Fig. 9(b) for c = 0.5, the (3 dependence remains weak, but the final 
transition probability |6(oo)p becomes very significant in general: it can even reach almost 
100% for certain values of (3. We hence observe nonlinearity-enhanced population transfer 
again, but now outside the RWA regime. As c is increased to c = 1.25 in Fig. 9(c), drastic 
and full-range oscillations of |6(oo)p are observed. Because Fig. 8(c) shows that this case 
corresponds to full Hamiltonian chaos, the full-range oscillation in Fig. 9(c) is clearly related 
to the underlying complete chaos on the mean-field level. As the nonlinearity parameter c 
further increases, the oscillation range of |6(oo)p decreases, consistent with the observation 
from Fig. 8 that the chaotic layer becomes smaller and smaller. In the case of c = 75 in Fig. 
9(f), there is essentially no transition probability for any value of [3. We have also done more 
detailed analysis of how mean-field trajectories should move as 7 is ramped, confirming that 
the size of the chaotic layer determines the oscillation range of |6(oo)|^. 
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FIG. 9: The final probability |6(oo)p as a function of phase /3 for a = 0.006 and to = 1. Results 
here should be connected with those in Fig. 8. 

Comparing results in Fig. 9(c) and 9(d) with those in Fig. 9(a) and 9(b), one finds 
another interesting aspect of driven NLZ dynamics. That is, contrary to our observation in 
the RWA regime, nonhnearity may also enhance the /9-dependence of |6(oo) p. Note, however, 
that the strong /3-dependence observed here is due to the mean-field chaos. Because chaos 
brings about strong instabilities and hence large fluctuations in the quantum system, results 
directly associated with the mean-fleld chaos should be applied with caution. For example, 
the validity of the mean-fleld treatment itself in the chaotic cases might be questionable when 
the number of particles in a BEC is not extremely large. This being so, a low-frequency 
driving fleld generating the mean-fleld chaos could be used to study effects of quantum 
fluctuations via the /9-dependence of |6(oo)p. 

VI. CONCLUSION 

Nonlinear extensions of LZ processes, and in particular, nonlinear LZ processes in a 
periodic driving fleld, are of importance in realizing the control of interacting many-body 
systems. By considering different parameter regimes of a periodic driving fleld, we have 
exposed a number of interesting dynamical phenomena in driven NLZ processes. Many of 
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these dynamical phenomena can be used for the control of NLZ processes by a driving field. 

In the high-frequency and large-amplitude regime, a driving field can effectively generate 
new self-interaction terms that are absent in non-driven two-mode NLZ models. The field 
can hence induce new energy-level structures and can be useful in simulating new systems 
not considered before. A driving-field-based scenario for the suppression of undesired NLZ 
tunneling probabilities also becomes possible. 

In the high-frequency and small-amplitude regime, driven NLZ processes are also much 
different from driven LZ processes in at least three aspects. First, for small nonlinearity 
strength the final transition probability shows a shifted phase dependence on the driving 
field. Second, for strong nonlinearity strength the transition probability shows a much 
weakened phase dependence upon the driving field. Third, in the adiabatic limit, significant 
transition probabilities become a purely nonlinear effect, because they were precisely zero 
in the absence of the self-interaction. These differences between driven NLZ and driven LZ 
processes suggest that a driving field may be used to directly expose and measure nonlinear 
effects. Furthermore, the adiabatic population transfer based entirely on the nonlinear 
nature of the system can be almost complete for certain amplitudes of the driving field, thus 
leading to an interesting approach to the improvement of adiabatic population transfer in 
the presence of nonlinear Landau-Zener tunneling. 

In the third low-frequency regime, nonlinearity may also enhance the final transition 
probability as well as its dependence on the phase of the driving field. The enhanced phase 
dependence occurs for intermediate nonlinear strength, where the NLZ dynamics displays 
Hamiltonian chaos. A driving field in this regime can therefore be used for studies of quantum 
fluctuations beyond the mean-field level, as well as studies of manifestations of mean-field 
chaos in the final transition probabilities. 

To conclude, we have exposed many interesting aspects of driven nonlinear Landau-Zener 
processes, by considering different parameter regimes of a single-frequency and constant- 
amplitude driving field. It is plausible that the driven nonlinear Landau-Zener processes 
are even much richer in the presence of more complex driving fields, such as those with 
multiple frequency components that have time-evolving amplitudes and phases. It should 
be also interesting to extend our nonlinear Landau-Zener treatments to other many-b ody 



generalizations of Landau-Zener processes that can go beyond the mean-field level 29 
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